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Simple  Algorithms  For  Calculation  Of  Scale  Factors 
For  Plane  Coordinate  Systems 
(1927  NAD  and  1983  NAD)  • 

by 


Dr.  Herbert  W.  Stoughton,  P.E.,  L.S. 


INTRODUCTION 

In  March  1972,  Berry  presented  a  paper  by  the  above 
title  at  the  Annual  Meeting  of  ACSM.  The  paper  discussed 
algorithms  for  the  state  plane  coordinate  systems  for 
1927  NAD.  Various  plane  coordinate  systems  (state  coordi¬ 
nate  systems,  city  systems.  Corps  of  Engineers  project 
systems,  the  Great  Lakes  system,  special  engineering  pro¬ 
jects,  etc.)  have  been  designed  for  the  purpose  of  simplify¬ 
ing  survey  calculations  from  the  spheroidal  earth  to  a 
plane  with  an  orthogonal  grid.  The  controlling  design 
criteria  in  many  instances  has  been  to  restrict  the  maximum 
difference  between  the  geodetic  length  and  the  grid  length 
to  one  part  in  ten  thousand  throughout  the  area  included 
within  the  system.  This  was  done  on  the  assumption  that  if 
the  reduction  of  measured  distances  to  geodetic  lengths  to 
grid  lengths  was  omitted,  the  surveys  would  meet  Third-Order 
accuracy  (one  part  in  five  thousand)  .  This  assumption  was 
valid  when  all  the  distances  were  measured  by  taping,  which 
was  normally  performed  to  0.01  or  0.02  feet  per  one  hundred 
feet,  and  thus  the  scale  factor  could  be  ignored. 

At  the  time  these  systems  were  designed,  this  assumption 
was  probably  valid,  but  supplementary  tables  were  required  to 
obtain  the  scale  factors  for  higher  order  surveys.  This 
meant  that  a  scale  factor  was  extracted  from  a  table  as  a 
function  of  latitude  (for  Lambert  systems)  or  the  perpendic¬ 
ular  distance  from  the  central  meridian  (transverse  Mercator 
systems).  With  the  advent  of  the  hand-held  computer/calcu¬ 
lator  and  digital  computers,  then  it  is  not  practical  or 
efficient  to  store  this  data  in  the  computer.  Furthermore, 
the  values  listed  in  these  tables  are  calculated  to  seven 
decimal  places,  which  is  insufficient  precision  for 
electronic  distance  measurements.  Also,  for  longer  lines 
it  is  recommended  that  the  scale  factor  for  a  line  should 
be : 


ri'J  XAB 
;  tr.urnouacfcd 
i  Jc-'tii‘ioati-n 

s» 

id  ! 

i 

i 

>> 

- , 

Pi3trit’u+ j.oa/ 

I  /vail  c: '  !  '  *  • -r  r  '■1?  3 


1 

A.  .Li  a 

.:.d/or 

i 

Diet. 

G^QC’ 

lal 

k  =  —  (k  +4 k  +  k  ) 

/\ 

6  i  m  2J 

H 

Where : 

=  the  scale  factor  at  one  end  of  the  line  (point  1). 

k  =  the  scale  factor  at  the  other  end  of  the  line  (point  2) . 

km  =  the  scale  factor  at  the  midpoint  of  the  line. 


Part  I 

[Iand-Held  Calculators  Version 


The  formulas  to  calculate  the  scale  factors  require  the 
geodetic  position  (latitude  and  longitude)  for  each  point. 

If  the  surveyor  is  employing  plane  coordinates,  this  requires 
a  transformation  of  the  plane  coordinates  to  geodetic  posi¬ 
tions  before  calculating  the  scale  factors.  The  formulas 
employed  to  calculate  scale  factors  appear  in  Part  II  of 
this  paper,  and  will  be  employed  to  derive  the  algorithms 
for  use  on  digital  computers. 


Derivation  Of  The  Algorithm 


Figure  1  is  a  section  depicting  the  relationship  between 
a  spherical  earth  and  the  projection  plane.  This  section  is 
taken  perpendicular  to  the  principal  axis  of  the  coordinate 
system,  which  for  various  projections  is: 

Lambert:  The  parallel  of  latitude  approximately 
midway  between  the  standard  parallels 

(designated  <p  )  ,  where  the  scale  factor 

o 

attains  its  minimum  value.  Then,  Figure  1 
is  a  meridional  section. 

Transverse  Mercator:  The  central  meridian  of  the 
zone  is  a  "normal"  section,  i.e.  it  is 
taken  perpendicular  to  the  central  meridian 
at  the  "approximate  center"  of  the  north- 
south  range  of  the  zone  [see  the  Appendix], 

and  designated  <f>  .  (This  does  not  trace 

o 

out  a  parallel  of  latitude.) 

Hotine  Skew  Orthomorphic  (Oblique  Mercator):  An 
approximate  geodesic  line,  which  passes 
through  a  specified  point,  and  has  a 
specified  azimuth  (non-cardinal  at  that 
point).  Figure  lisa  section  perpendicular 
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Figure  1 


to  this  line,  preferably  through  the 
specified  point  ("center  point"),  designated 


In  all  three  cases  the  spherical  radius  in  this  section 
is  designated  p  ,  and  is  calculated  by: 


a  J 1.  -  e2~ 

(1.  -  c2  sin2!?  ) 
0 


(1) 
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Where : 

a  =  the  length  of  the  semi-major  axis  of  the  reference 
ellipsoid  (in  meters) . 

ez  =  the  square  of  the  eccentricity  of  the  reference 
ellipsoid . 

Table  1  contains  the  numerical  values  of  a  and  e2  for 
1927  NAD  and  1983  NAD. 


Table  1 

Ellipsoid  a  (meters)  £2 

1927  NAD  6,378,206.400  0.006768657997 

1983  NAD  6,378,137.000  0.006694380025 


The  scale  factor  of  the  projection  along  the  principal 

axis  is  designated  k  .  At  any  point  P  located  at  a  distance 

o  * 

s,  perpendicular  to  the  projection's  principal  axis,  the 

scale  factor  is  k.  Disregarding  conformality  and  assuming 

spherical  shape  of  the  earth,  the  distance  from  0  (the  center 

of  the  earth)  to  point  Q  (the  point  on  the  projected  principal 

axis  where  it  is  intersected  by  the  depicted  section)  is 

eoual  to  the  radius,  o  ,  multinlied  by  the  central  meridian 

o 

scale  factor,  k  .  At  point  P,  the  distance  from  the  center 

o 

of  the  earth,  0,  is  p  k.  Then,  employing  the  Pythagorean 

o 

relationship  to  the  right  triangle  QPO  yields : 

(p  k)  2  =  (p  k  ) 2  +  s2 

0  0  0 

And : 

P  k  =  [(p  k  ) 2  +  sz  ]  ^  (2) 

0  *-oo  J 

Recall  the  binomial  series : 


,  ,  >.n  n  ,  n-1  .  n(n  -  1)  n-2  2 

(x  +  y )  =  x  +  nx  y  +  — i - -  x  y 

2 


,  n(n  -  1)  (n  -  2)  n-3  3  . 

+  — i _ L2 _ :  x  y  + 


=  (p  k  )2 
0  0 


Let : 


x 
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y  -  s2 

1 

n  =  T 


Then,  solving  for  k ,  equation  (2)  can  be  written: 


(P  )  + 

W  0 


—  (?/-:  )-:c2  -  — (p  a  )~2s4 

2  c  c  8  0  0 


+  —  (p  '<  ) 
16  0  0 


-5.5 


Rearranging : 


k  =  k 


Finally: 


2p2  a 

c  : 


S  P  L  •:  3 
o  o 


1 6  p £  a  5 
c  o 


K  = 


k  !  1 .  + 

n  I 


'  L 


i  2  p 2  k  2 


2  p  ‘  a  ^  /  2 

c  o 


I  1  i 


Up 


k'j  2  "’J 


0  5 


Let  : 


r2  =  io  -15  s 2 


A  = 


10 


1  5 


2p2*2 
0  0 


Where,  A  is  a  constant  for  each  system. 
Then,  equation  (3)  is: 


(3) 


k  =  k o  [l.  +  AT2(1.  -  -^r2(l.  -  Ar2  ...))]  (4) 

In  practice,  the  fourth  term  is  insignificant  and  can  be 
disregarded.  The  third  term  is  required  only  when  the  scale 
factor,  k ,  is  required  at  the  extreme,  width  of  the  system 
(greater  than  seventy  miles  from  the  principal  axis) . 

It  remains  to  develop  formulae  for  the  computation  of 
a ,  thus  T  (the  distance  from  the  principal  axis)  for  any 
point  P  defined  by  its  plane  coordinates  (x  -  east;  y  -  north). 
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Transverse  Mercator  Systems. 

r2  =  (x  -  x  )2  -no-1 5 

0 


(5) 


Where : 

x  =  the  x-coordinace  of  the  system's  central  meridian  (listed 
0  in  Table  2). 


Lambert  Conic  Systems. 


]:*10"15 


(6) 


'.There : 

a,  =  the 
o 

?.  =  the 


constant  for  each  system  [listed  in  Table  2], 
constant  for  each  system  [listed  in  Table  2). 


Kotine  Skew  Orthomorphic  (Oblique  Mercator)  Systems. 

T2  =  [(x  -  x  )  c  o  s  y  -  (u  -  ;/  )  sin  y  ]2  *10'15  (7) 

0  0  •  ‘  0  0 


Vine  re  : 

x  =  x- trans lation  constant  for  each  system  [listed  in  Table 

0  3]  . 

y  =  y- translation  constant  for  each  system  [listed  in  Table 
0  3], 

y  =  azimuth  of  the  principal  axis  at  the  center  point 
0  [listed  in  Table  3]. 


The  data  in  Table  2  refers  to  the  North  American  Datum 
of  1927.  A  similar  set  of  numerical  values  of  K  must  be 
calculated  for  1983  NAD.  At  this  time  (1981)  it  is  inappro¬ 
priate  to  prepare  similar  tables,  as  it  is  unknown  whether 
the  states  will  retain  the  current  systems  and  defining  par¬ 
ameters.  For  the  transverse  Mercator  systems  these  parameters 
(other  than  the  ellipsoid  and  datum)  the  central  meridian, 

X  ,  and  the  scale  factor  along  the  central  meridian,  k  ,  and 
0  o 

the  x-coordinate  of  the  central  meridian,  x  ,  need  to  be 

o 

defined.  Historically,  0  was  arbitrarily  selected  at  the 

o 

rnidgeodetic  latitude  of  the  zone,  rounded  to  the  nearest  five 

minutes  of  arc  [except  for  Alaska  and  Delaware].  If  the 

basic  design  of  the  existing,  transverse  Mercator  system  is 

identical  to  the  present  design,  £  listed  in  Table  2  is 

o 

acceptable  for  calculating  K.  Then,  for  1983  MAD  the 
ecju.it  ion  to  calculate  X  is: 


Constants  For  Lambert  And  Transverse  Mercator  Systems  (1927  NAD) 
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In  meters: 


12.3736955  (1.  -  0.006694380025  sin2!  ) 


k2 


(3-a) 


;  e  e  t : 


1-14955853  (1.  -  0.005694380025  sin2!  )2 

o 


(3-b) 


:he  numerical  values  of  ;;  in  Table  2  are  for  coordinates 


The  procedure  to  calculate  K  for  Lambert  systems  is 
identical,  but  simpler.  One  of  the  zone  constants  is  l . 
which  ('by  definition)  is: 


l  -  sin 


Then  ,  equations  (8-a)  and  (o-b)  are 


in  meters: 


„  _  12.3736955  (1.  -  0.006694380025  l2)2 

A  - -  ..  —  .  ■■■■■■■  — — -  — .  ■  ■  . .  ....  ■ — 

*2 

o 

In  feet: 


(8-c) 


1.14955853  (1.  -  0,006694380025  l2)2 

k2 


Accuracy  Of  Algorithm 

The  question  to  be  answered  is:  "How  good  is  the 
algorithm?"  A  comparison  was  made  between  the  scale  factors 
published  in  the  "projection  tables"  and  computed  by  the 
algorithm.  Table  4  is  a  comparison  of  the  scale  factors 
calculated  by  equation  (4)  and  published  in  the  projection 
tables.  California,  Zone  1,  was  selected  for  the  comparison. 
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Table  4 

Comparison  Of  Scale  Factors  -  Lambert  System 


(California,  Zone 

1) 

r 

k 

k 

differ. 

ratio 

r 

(£cet) 

(algorithn) 

(published) 

(*10~7) 

(l./diff .) 

(-.iles) 

0 

0.  S99S946 

0.9998946 

0 

- 

0.0 

6^6 

0. 999S946 

0. 999S946 

0 

- 

0.1 

31,004 

0. 5298557 

0.9356957 

0 

- 

5.9 

61,301 

0.  999S3S9 

0.99939SS 

0 

- 

11.6 

91,713 

0.9599042 

0.9995042 

1 

1:10,000,000 

17.4 

122,075 

0.9599117 

0.9999116 

1 

1:10,000,000 

22.1 

1 5  2  .  - :  3 1 

0.9599212 

0.9999211 

1 

1:10.000,000 

28.9 

1  S3. 7c 3 

0.9959328 

0.9999327 

1 

1:10,000,000 

3e  ,  6 

213,144 

0.9999466 

0.9999464 

2 

1:5,000,000 

40.4 

243,500 

0.9999624 

0.9999622 

2 

1:5,000,000 

46.1 

273,557 

0. 9999S04 

0.9999801 

3 

1:3,333,000 

51.9 

2C- , 213 

1.0300004 

1.0000000 

4 

1:2,500,000 

57.6 

334,570 

1.0000226 

1.0000220 

6 

1:1,667,000 

63.4 

364,927 

x. 0000469 

1.0000461 

8 

1:1,250,000 

69.1 

395,284 

1.0000732 

1.0000723 

9 

1:1,111,000 

74.9 

425,642 

1.0001017 

1.0001005 

12 

1:833,000 

80.6 

456,000 

1.0001323 

1.0001308 

15 

1:667,000 

86.4 

486,359 

1.0001650 

1.0001632 

18 

1:555,000 

92.1 

516,718 

1.0001998 

1.0001977 

21 

1:476,000 

97.9 

547,078 

1 .0002367 

1.0002341 

26 

1:384,000 

103.6 

The 

version  of 

equation 

(4)  employed  was: 

k 

=  k  (]  1 .  + 

o  *- 

AT2  (1.  -  0 

.5  AT2)] 

Where : 

k  =  0 . 9S98946353 
K  =  1.143216 
x  =  0 . 


y  =  y '  (from  projection  tables) 

At  the  extreme  limits  of  the  values  in  the  projection  tables 
(over  one  hundred  miles  from  the  principal  axis)  the  ratio 
of  the  difference  in  the  two  scale  factors  exceeded  one  part 
in  380,000.  At  the  geographical  limits  of  the  zone,  the 
ratio  of  the  difference  in  the  scale  factors  is  one  part  in 
600,000.  Incorporation  of  additional  (higher  order)  terms  in 
equation  (4)  did  not  alter  the  values  listed  in  Table  4. 

New  Mexico,  East  Zone,  was  selected,  and  equation  (4-a) 
was  employed.  For  the  range  of  values  listed  in  the  projection 
tables  the  scale  factor  was  calculated,  and  the  results  listed 
in  Table  5.  There  is  no  difference  between  the  numerical 
values  of  the  scale  factors.  Jhe  observed  difference  in  the 
seventh  decimal  place  can  be  attributed  to  numerical  roundoff, 
and  not  to  any  systematic  factors. . 


Table  5 


Comparison  Of  Scale  Factors  -  Transverse  Mercator  System 

(New  Mexico,  East  Zone) 


r 

k 

k 

differ'. 

7 

(feet) 

(algorithm) 

(published) 

(*10'7) 

(miles) 

0 

0.9999091 

0.9999091 

0 

0.0 

25,000 

0. 999909S 

0. 999909S 

0 

4.7 

50,000 

0.9999120 

0.9999121 

1 

9.5 

75,000 

0.9999155 

0.9999156 

1 

14.2 

100,000 

0.9999205 

0.9999205 

0 

18.9 

125,000 

0.9999270 

0.9999270 

0 

23.7 

150,000 

0.9999343 

0.9999348 

0 

2S.4 

175,000 

0.9999441 

0.9999441 

0 

33.1 

200,000 

0.9999549 

0.9999549 

0 

37.9 

225,000 

0.9999670 

0.9999671 

1 

42.6 

250,000 

0.9999-306 

0. 9S99S06 

0 

47.3 

275,000 

0.9999957 

0.9999956 

1 

52.1 

300,000 

1.0000121 

1.0000121 

0 

56.8 

325,000 

1.0000300 

1.0000300 

0 

61.6 

350,000 

1.0000493 

1.0000493 

0 

66.3 

375,000 

1.0000701 

1.0000701 

0 

71.0 

400,000 

1 .0000922 

1 .0000922 

0 

75.8 

425,000 

1 .  Ouu  1 ' 

1.0C01 15S 

0 

SO.  5 

450 , 000 

l.OOuliOo 

1.0001409 

1 

85.2 

475,000 

1.0001673 

1.0001674 

1 

90.0 

500,000 

1.0001952 

1.0001953 

1 

94.7 

525,000 

1.0002246 

1.0002246 

0 

99.4 

530,000 

1.0:02306 

1  .0002  397 

0 

100.4 
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Application  To  UTM 


In  the  previous  discussion,  the  algorithm  was  applied 

to  systems  with  relatively  narrow  north-couth  extent. 

Excluding  Alaska,  the  maximum  change  in  latitude  is  about 

three  degrees  of  arc.  The  UTM  systems  extend  over  the  range 

oFHP#J)  to  80°S.  Therefore,  a  single  numerical  value  of  K 

may  net  be  appropriate.  Table  6  lists  the  numerical  values 

of  X  for  every  2.5  degrees  of  latitude  for  the  existing 

sysLems  for  metric  coordinates  [Clarke  1866  ellipsoid  and 

k  =  0.9996000000} . 

0 


Tab!  e  6 

Values  Of  K  For  Intervals  Of  2.5  Degrees  Of  Latitude 


(Clarke  1866  Ellipsoid) 


[metric 

coot  di nates ) 

K 

c 

K 

K 

K 

0.0 

12.3342570 

32.5 

12.3359054 

65.0 

12.2469335 

2.5 

12 . 3S393SO 

35.0 

12.3291634 

67.5 

12.2415725 

5.0 

12.3323535 

T  -  r 

12.3122056 

70.0 

12.2366611 

7.5 

12 . 3816009 

40.0 

12 .3150651 

72.5 

12.2322363 

10.0 

12.3792 023 

42.5 

12 . 3078562 

75.0 

12.2283317 

12.5 

12.3764045 

45 . 0 

12.3005740 

77.5 

12.2249766 

15.0 

12.3730291 

47.5 

12.2932940 

80.0 

12.2221963 

17.5 

12.3691021 

50.0 

12.2860716 

82.5 

12.2200119 

20.0 

12.3646535 

52.5 

12.2789615 

85.0 

12.2184397 

22.5 

12.3597174 

55.0 

12.2720179 

87.5 

12.2174916 

25.0 

12.3543318 

57.5 

12.2652934 

90.0 

12.2171748 

27.5 

12.3485379 

60.0 

12.2588389 

30.0 

12.3423801 

62.5 

12.2527035 

The  value  of  } 

'  changes 

with  latitude 

in  a 

nonlinear 

manner,  but  is  relatively  insensitive  to  small  changes  (a 
few  tens  of  minutes  of  arc)  in  latitude. 


Let: 


A.nd : 


10 


1  5 


6  2o2(l.  -  e2)k2 


(9) 


0  =  _ _ .  (in  radians) 

P* 


(10) 
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Where  : 


175  C«1 
2  56  J 

(ID 


i  .  _  '  . 

K  K  ( 1 .  -  eisin2C)1  (12) 

0 

T.~  .  r  . 

a  the  se  mnrs  a  jor  sxis  cf  the  ellipsoid  (in  peters). 

g1  the  squaie  of  f he  eccentricity  of  the  ellipsoid. 

k  the  scale  factor  along  the  central  meridian  of  the 
0  projection  system 

for  the-  Clarke  186  6  ellipsoid  (1  327  NAD),  equation  (12) 
[for  metric  coordinates)  (k  -  0  9996CGCC0C)  is-. 


9"  -  0 . 03  2  3  9.38  j  JL 

k 

c 

[ihte  s  *-"  is  the  r.ngl  •  Q  in  second  ■  of  ana.] 

K  12*84257  (  1.  -  C. 006768657997  sin3e)2 

(12-a) 

For  the  1983  i.'AD  ellipsoid,  equation  (12)  [for  metric 

coordinates]  (<  =  0 . 999o0000GG)  is: 

o 

0"  =  0.03239363  JL 

k 

0 

X  -  12.383600  (1.  -  0 . 005 696 380025  sin2©) 2 

(12-b) 

If  the  project  does  not  extend  more  than  about  two 
degrees  in  latitude,  select  a  mean  y-coordinate  for  the  pro¬ 
ject  to  calculate  0.  There  will  be  no  Iocs  in  fne.jirion. 

Furthermore,  to  ret  3 in  a  precision  of  one  in  the  ninth 
decimal  place  of  the  computed  scale  factor,  the  last  term 
given  in  equation  (4)  [-  Tf T ~  ]  must  be  incorporated  for  dis¬ 
tances  greater  than  250  kilometer.,  from  the  central  meridian. 
For  shorter  distances  from  t:'.n-  central  meridian,  equation 
(4-a)  is  sufficient. 
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Fart  2 

Digital  Computer  Version 


The  previous  material  discussed  an  algorithm  for  com¬ 
puter/calculator  with  limited  data  storage  capability.  With 
programmable  computers,  an  algorithm  can  be  developed  to 
employ  the  plane  coordinates  and  the  published  equations  that 
require  the  geodetic  positions.  These  algorithms  will  con¬ 
vert  the  given  coordinates  to  requisite  geodetic  values  for 
computing  the  scale  factors. 


Transverse  Mercator  Systems 


The  formula  to  calculate  scale  factors  from  a  given  set 
of  plane  coordinates,  x  and  y ,  requires  the  footpoint  geodetic 
latitude,  <J>^.  The  footpoint  latitude  at  the  point  is  calcu¬ 
lated  from  the  y-coordinate .  To  calculate  <j>  requires  the 

meridional  arc  length  from  the  equator.  Except  for  UTM 
systems,  the  origin  of  the  ^-coordinates  ( y  =  0)  is  a 
parallel  of  latitude  between  the  equator  and  the  "working" 

area  of  the  system.  This  requires  a  constant,  y  ,  for  each 

o 

zone.  To  obtain  tne  numerical  value  of  y  (for  the  precision 
required),  let: 


A  =  1.  +  JL  e2 

+  45  e-  +  175  e* 

0  4 

64  256 

A  =  1  +11 

c-  +  125  es 

1  4  16 

512 

A  =  15  4- 

105  es 

2  128 

512 

A  =  — e  6 

3  1536 

Then  (in  meters) : 

H  ~  a (1 .  -  c 2 )  [  A  )  +  A  sin  2$  +  A  sin  4J 

o  l  o  r  i  2 

+  A^sin  6^1 


(13) 
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Where : 

a  -  the  semimajor  axis  of  the  ellipsoid  (in  meters). 
e2  =  the  square  of  the  eccentricity  of  the  ellipsoid. 

$  =  the  geodetic  latitude  of  the  origin  of  the  coordinate 

system.  Note:  the  subscript  r  designates  the  angle  in 
radians . 

The  value  of  y  is  a  constant  for  a  zone.  The  next  step 
o 

is  to  calculate  the  footpoint  latitude.  This  is  performed  in 
two  steps. 

e"  =  v(~y  +  yj  (i;> 

k 

0 

Where:  y  and  y  are  in  meters,  and  y  is  a  constant  for  an 
o 

ellipsoid  (listed  in  Table  7) .  0"  is  the  approximate  value 
of  the  latitude  in  seconds  of  arc.  Then,  the  footpoint 
geodetic  latitude,  ,  is: 

*  -  0  +  3  sin  20  (15) 

Table  7 

Ellipsoid  y  3 

Clarke  1866  (1927  NAD)  0.03239388  52573 

1983  NAD  0.03239363  51975 

The  units  of  the  term  "8  sin  26"  in  equation  (15)  is  seconds 
of  arc. 

Then,  let: 


t  =  tan  y 

i  i 


K  *  .* 

C  1 


Then,  the  scale  factor  is: 


*  •  *■  I1'  +  I  f 


1.  +  —  |(l.  +  n2)  +  —  (1.  +  n 2  ( 6 .  + 

2  1  12  1 


9nJ  +  4r, ’  -  24t2ry(l.  +  q 

i  i  ’  l 


*»r 

i 


The  precision  of  equation  (16)  is  sufficient 
for  the  UTM  systems  as  well  as  the  state  plane  coordinate 


lA.".BER  i  SySTz'-'S 

The  procedure  for  calculating  the  scale  factor  utilizes 
a  different  set  of  formulas,  but  employs  an  analogous  pro¬ 
cedure.  For  each  Lambert  system,  the  following  constants  are 

provided:  x  ,  l  (sin  p  )  ,  and  X .  Note:  the  X  just 

D  0  0 

mentioned  is  not  the  constant,  X,  that  was  employed  earlier 
to  calculate  the  scale  factors  in  Part  1  of  this  paper. 

Then  : 


6  =  arc tan 


p  -  2  arc  tan 


H.  -  u 

o 


i  s  in  G 


10  eg 2 

9 . r 2o 2736 


1.  + 


•:  lit . 


The  units  of  the  second  term  in  this  equation  is  s 
c Then,  the  scale  factor  is; 


l  X  71.  -  a  2  s  i  n 2  o 


The  tend  ;tic  latitude,  employed  in  equal. ion  (17)  will  be 
ac.  • . r . n  ..  c o  bettor  t!;tn  ;  ive  secon  it  of  .arc.  This  is 
.n.iT  •  c :  :.t  lot  .  :>•  anting  >  to  air.,  r  _c::c  il  r!  . 
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Appendix 


Calculation  Of  Jj  For  The  State  Plane 

o 

Coordinates  -  Transverse  Mercator  Systems 


by 

Dr.  Herbert  W.  Stoughton,  P.E.,  L.S. 


^  *-ri  icicle  Or  Constants"  ieLrber  oe^e  7  c**  10 

cr  the  "projection  tables"  for  the  state]* is  listed*  toe 

'll..  '  ' . . 


nvrr.art; 


al  ur- 1  "1 


i0o  I 


L  6  p  2  j 


mxs  number  is  a  function 


raus  to  iccat 


me 


of  tne  geodetic  latitude  i  . 

Wines  ot  ror  all  the  transverse  Mercator  systems  in  the 

Lnited  States.  This  Appendix  contains  a  derivation  to 
calculate  based  upon  the  numerical  values  listed  in  the 

aforementioned  "Table  Of  Constants". 

The  symbol  p#  is  the  "mean  radius"  of  the  ellipsoid 
for  geodetic  latitude  <J>  . 


Where : 


p.  ■  AvT 
0 


(A- 1) 


R 


«(1.  ~  ^2) 


(I-  -  c2.mn2:  ) 


3/2 


(A- 2) 


11  = 


•  -  a  2  sin2  ) 


1/2 


(A- 3) 


C  =  fRuatorial  radius  (wvimajor  axis)  cf  the  ollinscid  (in 
meters).  •  v  " 

e2  -  the  square  of  the  eccentricity  of  the  ellipsoid. 

Substituting  equation  (A- 2)  and  (A- 3)  into  ecuaticn  (A-I 


0  ■"  (  '  . 


) 


■  _  1  -- 


:1:1  V 


a. 
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The  formula  to  calculate  "log 


'  1  „ 

6  o2 
•  « 


[based  upon  U.S.C.&G.S. 


Special  Publication  Nos.  193  and  195]  is: 


r  i  l 

os  !  - —  !  = 


los  i  rr' 

l6pJ 


(1.  -  a  2  sin2  q  )2  / 1200 1 2 


k2  e2(l.  -  e2) 
o 


3937 


-r  20. 


(A-4) 


=  the  scale  factor  along  the  central  meridian  of  the  system 


The  constant  20.  is  added  to  obtain  a  positive  characteristic. 
For  the  Clarke  lo65  ellipsoid  (1927  NAD) : 

a  =  6,378,206.400  meters 
e 2  =  0.0067686580 
Then,  equation  (A-4)  is: 

log  JL  =  log  f  3 . 83206*10  ~ 1 5  (1-  -  a  2  sin2  +  2o. 

.  6p2  .  fc2 


Let: 


l[»  =  lo 


^  r  1 

.  6p‘ 


(A- 4a) 


equation  (A- 4a)  can 


rp  =  log  3. 8  3206  “10  : 3  4-2  log  (1.  -  c2sin2y  ) 

-  2  log  4-  20.  (..-5) 

c 

Th  2  numeric::!  wit:  a  o'  y  i r  Ivan  in  the  "Tab  la  Of 
•ir.ts".  Rearranging  equation  (A- 5)  yields: 


!-:■  (1. 


2  sin2  :  )  -  —  :  -  leg  3 .  o32C-6*'flC*‘ ‘ 

0  O  ; 


+  2  lo:: 


(-.  ; 
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Equation  (A-6)  contains  the  known  quantities  [ ^  and 

on  the  right  side  of  the  equation,  and  the  unknown  quantity, 
i  .  is  cn  the  left  side  of  the  equation. 

La  t : 


leg  a  -  log  3 .  S3206'-i0  “IS  +  2  log 

2  i 

Then,  equation  (A-6)  can  be  written: 


20 


log  (1.  -  a ;  sir.  S'.) 


(A- 7) 


-.r.  c : 


2  ■  *>  x 

e  sm~ 6 

o 


V 


Finally : 


(A- 8) 
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